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ABSTRACT 

Massive Black Hole (MBH) seeds at redshift z > 10 are now thought to be key ingredients 
to explain the presence of the super-massive (10 ^“**'Mq) black holes in place < IGyr af¬ 
ter the Big Bang. Once formed, massive seeds grow and emit copious amounts of radiation 
by accreting the left-over halo gas; their spectrum can then provide crucial information on 
their evolution. By combining radiation-hydrodynamic and spectral synthesis codes, we sim¬ 
ulate the time-evolving spectrum emerging from the host halo of a MBH seed with initial 
mass 10^ Mo, assuming both standard Eddington-limited accretion, or slim accretion disks, 
appropriate for super-Eddington flows. The emission occurs predominantly in the observed 
infrared-submm (1 - 1000jum) and X-ray (0.1 - 100keV) bands. Such signal should be easily 
detectable by JWST around ~ 1 jim up to z ~ 25, and by ATHENA (between 0.1 and lOkeV, 
up to z ~ 15). Ultra-deep X-ray surveys like the Chandra Deep Eield South could have already 
detected these systems up to z ~ 15. Based on this, we provide an upper limit for the z > 6 MBH 
mass density of p, <2.5 x 10 ^MqMpc“^ assuming standard Eddington-limited accretion. If 
accretion occurs in the slim disk mode the limits are much weaker, p, <7.6 x lO^MoMpc"^ 
in the most constraining case. 

Key words: accretion - black hole physics - quasars: supermassive black holes - radiative 
transfer - cosmology: dark ages, reionization, first stars - cosmology: early Universe 


1 INTRODUCTION 

The cosmic epoch in the redshift range 10 < z < 30 was charac¬ 
terized by the formation of the first stars (Pop III) and of the first 
black holes (see Bromm & Yoshida 2011, Volonteri 2010, Volon- 
teri & Bellovary 2012, Haiman 2013 for recent reviews). Detecting 
these sources directly is at, or beyond, the sensitivity edge of cur¬ 
rent observatories. While recently Sobral et al. (2015) have shown 
that a very luminous Lyo- emitter at z x: 6.6 may be consistent with 
having a mixed composition of Pop III and second-generation (Pop 
II) stars, to date there are no confirmed observations of the first 
black holes, partly due to the uncertainty on their observational 
signatures. The next generation of observatories will most likely 
detect the first glimpses of light in the Universe, both in the elec¬ 
tromagnetic spectrum (e.g. ALMA, JWST, ATHENA) and in the 
gravitational waves domain (e.g. e-LISA, DECIGO). 

The formation process of the first black holes is likely to pro¬ 
duce a strong imprint on their observational signatures, as well as 
on their mass growth. The standard theory of Eddington-limited 
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accretion predicts that black holes grow in mass over a time scale 
~ 0.045eo.i Gyr, where eo.i is the matter-energy conversion factor 
normalized to the standard value of 10%. With black hole seeds 
of initial mass ~ IOOMq, formed at the end of the very short 
(~ 1 - lOMyr) lifetime of Pop III stars, it is, at best, challeng¬ 
ing to explain recent observations of optically bright quasars with 
M. ~ 10®^'°Mo at z ~ 7 (Mortlock et al. 2011; Wu et al. 2015). 
An alternative, attractive solution is based on Massive Black Hole 
(MBH) seeds (10^“^ Mq) appearing at z ^ 10- 15, giving a jump 
start to the growth process (e.g., Spaans & Silk 2006; Begelman 
et al. 2006; Lodato & Natarajan 2006). Under specific conditions 
(Bromm & Loeb 2003; Begelman et al. 2006; Volonteri et al. 2008; 
Shang et al. 2010; Johnson et al. 2012; Agarwal et al. 2014), the 
collapse of a primordial atomic-cooling halo may lead to the forma¬ 
tion of MBHs with a birth mass function peaked at Af, ~ 2x 10^ Mq 
(E errara et al. 2014). The subsequent gas accretion from the host 
halo leads to further growth into > 10^ Mq objects. 

The expected abundance of MBHs is largely unconstrained. 
Salvaterra et al. (2012) and Treister et al. (2013) provided upper 
limits of order p, < lO^^'^MoMpc^^, using the X-ray background 
and the stacked X-ray luminosity of high-redshift galaxies, respec¬ 
tively. Upper limits for the z = 6 MBH mass density provided by 
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Willott (2011), Fiore et al. (2012) and Cowie et al. (2012) are even 
lower (p, < lO^^^MoMpc'^), although none of these constraints 
take into account Compton-thick sources that can be buried deep 
inside dense nuclei in proto-galaxies. Yue et al. (2013, 2014) noted 
that if MBHs are responsible for the near-infrared background fluc¬ 
tuations, their high-redshift mass density should be comparable to 
the present-day value: p,(z = 0) ^ 2x 10^ MoMpc^^, see Yu & 
Tremaine (2002). The population of high-redshift MBHs would 
produce gravitational waves, during their collapse (Pacucci et al. 
2015b) or ensuing MBH-MBH mergers (Sesana et al. 2007, 201 1), 
detectable with upcoming observatories. 

Our previous works (Pacucci & Ferrara 2015; Pacucci et al. 
2015a) focused, through accurate ID radiation-hydrodynamic sim¬ 
ulations, on the dynamical evolution of z ~ 10 MBH seeds with 
initial mass M, ~ 10^“®Mq, embedded in dark matter halos with 
total mass M/, ~ 10^ Mq and accreting in the standard Eddington- 
limited scenario, or including a model for super-Eddington accre¬ 
tion through slim disks (see also Volonteri & Rees 2005; Volon- 
teri et al. 2015). In the present work we focus on their emission 
spectrum with three objectives: (i) predict the time evolution of the 
spectrum, (ii) assess the observability with current (Chandra Deep 
Field South, CDF-S) and future (JWST, ATHENA) surveys, and 
(iii) estimate the mass density p, of high-redshift MBHs. 

The outline of this paper is as follows. In §2 we de¬ 
scribe the physical and numerical implementation, while in 
§3 we present our results for the spectral evolution of high- 
redshift MBHs. Finally, in §4 we provide some further dis¬ 
cussion and the conclusions. We adopt recent Planck (Planck 
Collaboration et al. 2015) cosmological parameters throughout: 
(n„„DA,ni,l!,«i,o-8) = (0.32,0.68,0.05,0.67,0.96,0.83). 


2 PHYSICAL AND NUMERICAL IMPLEMENTATION 

The present work is based on radiation-hydrodynamic simulations 
post-processed with CLOUDY, a spectral synthesis code (Ferland 
etal. 2013). 

The physical framework is the following: a high-redshift (z = 
10) MBH seed with initial mass 10^ Mq is located at the center of a 
dark matter halo with primordial composition and total mass Af/, ~ 
10^ Mq (Tyir ~ lO'* K). The MBH accretes mass from the inner parts 
(within ~ lOpc) of the host halo until complete gas depletion. 

Our radiation-hydrodynamic code takes into account the 
frequency-integrated radiative transfer through the gas, with appro¬ 
priate: (i) cooling and heating terms, (ii) matter-to-radiation cou¬ 
pling, and (iii) energy propagation through a two-stream approxi¬ 
mation method. The code computes the accretion rate through the 
inner boundary of the simulation domain, from which we derive 
the total bolometric energy radiated by accretion, assuming a ra- 
diatively efficient or inefficient disk. The f\i\\ frequency-dependent 
radiative transfer through the host halo is then performed in a 
post-processing step using CLOUDY. This code computes the de¬ 
tailed time-evolving spectrum emerging from the host halo us¬ 
ing as input the matter distribution obtained from our radiation- 
hydrodynamics simulations and the realistic irradiation spectrum 
at the inner boundary, scaled to the appropriate bolometric lumi¬ 
nosity. Additional details are given in sections 2.1 and 2.2. 

2.1 Dynamics and thermodynamics 

Our radiation-hydrodynamic code (see Pacucci & Ferrara 2015 
for an extensive description) solves the ID spherically-symmetric 


equations of hydrodynamics and a frequency-integrated version of 
radiative transfer equations. The code evolves self-consistently the 
radial component of the standard system of ideal, non-relativistic 
Euler’s equations (neglecting viscosity, thermal conduction and 
magnetic fields) for a gas accreting, with no angular momentum, 
onto the central MBH, assumed at rest and already formed at the 
time t = Q, with a given initial mass M,{t = 0). The simulation do¬ 
main spans from 0.1 pc to lOpc, largely encompassing the charac¬ 
teristic spatial scale for accretion, the Bondi radius: 

~ 3.0pc, (1) 

c / X 

5 ( 00 ) 

where G is the gravitational constant and Cy(oo) = ^lyRToclp ~ 
12kms“* is the sound speed at large distances from the accretion 
boundary; y = 5/3 is the ratio of specific heats, R is the gas con¬ 
stant, T is the gas temperature and /r = 1.15 is the mean molecular 
weight for a primordial H-He composition gas with helium fraction 
Yp = 0.24665 (Planck Collaboration et al. 2015) and no metals. For 
a 10^ Mq object, the inner boundary of our spatial domain is ~ 10^ 
times larger than the Schwarzschild radius and ~ 10^ times larger 
than the centrifugal radius, i.e. the spatial scale below which de¬ 
viations from spherical symmetry become important and an accre¬ 
tion disk may form. Moreover, the angular momentum transfer in 
the outward direction of the accretion flow is very efficient, due to 
gravitational torques induced by dark matter and gas distributions 
of the halo (Choi et al. 2015). The gas loses its angular momentum 
efficiently and flows well beyond its centrifugal barrier. Therefore, 
despite its simplifications, our ID approach is significantly help¬ 
ful in acquiring physical insights on the mechanisms regulating the 
black hole growth. For the same reason, neglecting viscosity, ther¬ 
mal conduction and magnetic fields is a safe choice, since they play 
an important role only on spatial scales comparable with the radius 
of the accretion disk. 

The gas accretion through the inner boundary of our spatial 
domain produces an accretion rate M,, which in turn generates 
an emitted luminosity L via two different accretion models: (i) a 
standard Eddington-limited model in which L - ec^M, e = 0.1, 
and (ii) a radiatively inefficient model, the slim disk, in which 
L oc ln(M,) and e <0.04 is a function of M,. The main physical 
quantity which determines the properties of the accretion disk, and 
consequently the radiative efficiency, is the accretion rate. Accre¬ 
tion of gas at moderate rates (0.01 <fEdd < 1) is expected to form a 
radiatively efficient, geometrically thin and optically thick accre¬ 
tion disk, which is typically modelled with the standard a-disk 
model (Shakura & Sunyaev 1973). In a Shakura & Sunyaev disk 
the radiative efficiency is determined only by the location of the 
innermost stable circular orbit, which in turn depends only on the 
spin of the black hole, and it varies between ~ 6% and ~ 32% for 
a non-spinning and a maximally spinning black hole, respectively 
(Thorne 1974). In a super-critical (fpdd > 1) accretion environment 
the structure of the accretion disk is, instead, expected to be ge¬ 
ometrically and optically thick and radiatively inefficient (but see 
Jiang et al. 2014; McKinney et al. 2015). The most common solu¬ 
tion proposed for such accretion flows is the slim disk (Paczynski & 
Abramowicz 1982; Abramowicz et al. 1988; Mineshige et al. 2000; 
Sadowski 2009, 2011; McKinney et al. 2014). In these scenarios, 
part of the energy produced inside the disk is advected inwards 
(see, e.g., Abramowicz & Eragile 2013; Lasota 2015) out to a spa¬ 
tial scale named the photon trapping radius (R,r ~ Rsfsdd, where 
Rs is the Schwarzschild radius). Therefore, only a fraction of the 
photons produced in the accretion disk is able to free stream out 
of Rfr- consequently, the effective radiation and radiation pressure 
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escaping to infinity is decreased (see e.g. Begelman 1978; Ohsuga 
et al. 2002). While the slim disk solution is the simplest and most 
tested model for super-critical accretion, alternatives exist, e.g the 
ZEro-BeRnoulli Accretion (ZEBRA, Coughlin & Begelman 2014) 
and the ADiabatic Inflow-Outflow Solutions (ADIOS, Blandford & 
Begelman 1999; Begelman 2012) models. These theoretical models 
also include a parameter which describes the fraction of the inflow¬ 
ing mass which is lost due to radiation pressure. 

Radiation pressure accelerates the gas via: 




K(p,T)Ur) 

AtiEc 


( 2 ) 


where the gas opacity k(p,T) includes Thomson and bound-free 
terms, with the inclusion of a temperature dependence (Begelman 
et al. 2008). The radiation pressure may be able to temporarily 
interrupt the gas inflow, resulting in an intermittent accretion and 
outflows. The physical parameters regulating this occurrence are 
investigated in Pacucci et al. (2015a). 

We assume that the gas initially follows the isothermal (T ~ 
lO'^K) density profile derived from the simulations in Latif et al. 
(2013a), approximated by the functional form: 


p(r) = 


PO 

1 -I- {rlaf ’ 


(3) 


where a is the core radius and po is the central density. To under¬ 
stand how the matter distribution influences both the accretion and 
the emerging spectrum (through the hydrogen column density), we 
implemented two different density profiles, both of them yielding a 
gas mass ~ 10^ Mq: (i) a high density profile (HDP) with a central 
density po = 10^'^gcm^^ and a core radius a = 0.002pc, and (ii) a 
low density profile (LDP) with a central density po = 10^'^gcm^^ 
and a core radius a = 2pc. The LDP may be thought as the density 
profile resulting after the formation of a MBH of mass ~ 10^ Mq 
at the halo center (see Latif et al. 2013b, 2014), while in the HDP 
case the seed formed is very small (< 10^ Mq). A density floor of 
lO^^'^gcm^^ ~ lcm“^, a factor at least 10® times smaller than the 
central density for both profiles, is imposed for numerical stabil¬ 
ity reasons. In summary, we have four models: standard accretion 
- LDP, standard accretion - HDP, slim disk accretion - LDP, slim 
disk accretion - HDP. 


2.2 Spectrum 

Our code evolves the system in time until gas depletion and pro¬ 
vides CLOUDY with data to compute the spectrum emerging from 
the host halo. This code takes into account: (i) the spatial profiles 
for hydrogen number density nup) and temperature T(r), (ii) the 
source spectrum of the central object, and (iii) the bolometric lu¬ 
minosity L of the source, computed self-consistently from M,. The 
spherical cloud of gas is assumed to be metal-free, so that only H 
and He recombination lines are present. The addition of a small 
amount of metals, formed by the first Pop III stars, would increase 
the number of lines and the absorption of high-energy photons. The 
source spectral energy distribution is taken from Yue et al. (2013) 
and can be described as the sum of three components: (i) a multi¬ 
colour blackbody, (ii) a power-law, and (iii) a reflection component. 
The source spectrum is extended from far-infrared to hard X-ray 
(~ IMeV). 


3 RESULTS 


We describe in the following the standard accretion - LDP case, and 
discuss the differences with the other cases when needed. 

The balance between the inward gravitational acceleration and 
the outward radiation pressure keeps the accretion rate close to the 
Eddington level (fzdd ~ L2) for most of the time, ensuring a contin¬ 
uous accretion. The top panel of Eig. 1 shows the time evolution of 
fsdd and of 3/. in this accretion scenario. Defining the duty cycle, 
2), as the fraction of time spent accreting during the total evolu¬ 
tionary time of the system, we find that in the LDP cases 0=1, 
while other accretion scenarios may be characterized by quiescent 
phases (O < 1, see Table 1). The system evolves for ~ 120Myr, un¬ 
til complete depletion of the ~ 10^ Mq gas reservoir within ~ lOpc. 
However, this time scale is a lower limit, since in a real galaxy the 
gas would extend much further. 

Fig. 2 (top panel) shows the time evolution of the bolometric 
luminosity emitted at the inner boundary, before traversing the host 
halo gas. The corresponding Eddington luminosity: 


l-Edd = 


4nGM,c 
-cv M ,, 

Kj 


(4) 


with Kj being the Thomson opacity, is also shown for comparison. 
The luminosity increases for the first ~ 115Myr, reaching a peak 
of ~ 5 X lO^^ergs^* and is, on average, mildly super-Eddington 
as long as the amount of gas is sufficient to sustain this accretion 
rate. Afterwards, the luminosity plummets, when all the available 
gas has been consumed. For comparison, the interested reader is 
referred to Pacucci & Ferrara (2015), where a plot of the same 
physical quantities is shown in Fig. 8. Moreover, Fig. 2 also shows 
the values of the hydrogen column density, Nh, computed at se¬ 
lected times {ts = 5,75, 110, 115, 120Myr) when also the emer¬ 
gent spectrum is computed. The system is initially Compton-thick 
{Nh > 1.5 X lO^^cm^^, see the horizontal line). As the gas is pro¬ 
gressively accreted by the MBH, Nh steadily decreases before a 
sudden drop {ts ~ 115 Myr) corresponding to the remaining gas 
evacuation by radiation pressure. Table 1 lists duty cycles and ac¬ 
cretion time scales for all four models. 

The bottom panels of Fig. 1 and Fig. 2 show the time evolution 
of fEdd^ tVl,, L and Nh in the slim disk - LDP case. In this accre¬ 
tion scenario, the Eddington rate reaches high values {fEdd ~ 20) 
only for a short amount of time (^ 2 Myr), while afterwards the ac¬ 
cretion rates are sub-Eddington. The evolution is much more rapid 
(the available gas is consumed in ~ 6Myr, a factor of ~ 16 faster 
than in the standard case) and the host halo becomes Compton-thin 
{Nh < L5x lO^'^cm^^) in < 2Myr. As already noted in Pacucci 
et al. (2015a), in the slim disk accretion - LDP case, the black hole 
is able to accrete up to ~ 80% of the gas mass within ~ lOpc, with 
respect to the ~ 15% in the standard accretion scenario, and the 
black hole grows to ~ 8 x 10® Mq in ~ 6 Myr, while in the standard 
case the black hole grows to ~ 1.5 x 10® Mq in ~ 100 Myr. 


3.1 Spectral evolution 

The time-evolving spectrum emerging from the host halo, in the 
standard accretion - LDP case, is shown in Fig. 3 at z = 9, i.e. 
~ 100Myr after the beginning of the simulation {z = 10). The spec¬ 
trum is composed by: (i) the continuum emitted by the source and 
attenuated by the gas, and (ii) the diffuse emission of the gas. Most 
of the energy emerges in the observed infrared and X-ray bands. 
The latter is characterized by a bell-shaped spectrum peaked around 
IkeV, while in the infrared band a large number of H-He neb¬ 
ular lines is present. Photons with frequency shortwards than the 
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Figure 1. Time evolution of the Eddington factor and of the black hole 
mass M,, in the standard accretion - LDP (top) and the slim disk accretion 
- LDP (bottom) scenarios. The Eddington ratio is reported as a running 
average over periods of ~ 0.1 Myr. 


Figure 3. Time evolution of the spectrum emerging from the host halo for a 
source located at z = 9, in the standard accretion - LDP case. The infrared, 
optical and X-ray bands are highlighted with shaded regions, while the un¬ 
processed spectrum is reported, at peak luminosity (t = 115 Myr), with a 
dashed line. The flux limits for future (JWST, ATHENA) and current (CDF- 
S) surveys are also shown. 



E 



Figure 2. Time evolution of the holometric luminosity emitted at the inner 
boundary, before traversing the host halo gas, in the standard accretion - 
LDP (top) and the slim disk accretion - LDP (bottom) scenarios. The lu¬ 
minosity is reported as a running average over periods of ~ 0.1 Myr. The 
corresponding Eddington luminosity, oc M., and the values of the hy¬ 
drogen column density (blue curves, right axis) are also shown. 


Lyff line are absorbed by the intervening matter at column densi¬ 
ties Nh > lO^^cm^^ and reprocessed at lower energies, boosting 
the infrared emission of the halo. X-ray emission occurs predom¬ 
inantly within the rest-frame energy range 4keV < Ey < lOkeV. 
The mean free path of such photons is much larger than the Hub¬ 
ble radius at z = 9: Tx > 8Gpc » Rh{z = 9) = 165 Mpc. Hence 
growing MBH seeds negligibly contribute to reionization. The in¬ 
crease with time of the continuum normalization in the X-ray is 
mainly due to the progressive rise of the bolometric luminosity of 


the central object. The ratio between the infrared and the X-ray 
continua depends on the column density, since in the Compton- 
thick case {Nh > 1-5 x lO^^cm^^) the high-energy frequencies are 
heavily absorbed and reprocessed at lower energies, leading to an 
overall increase of the infrared emission. When the gas becomes 
Compton-thin at ~ 75 Myr, the X-ray continuum progressively in¬ 
creases, while the infrared one starts to decrease. Approaching the 
complete gas depletion within ~ lOpc, at ~ 120Myr, the column 
density is so low (see Fig. 2, top panel) that the outgoing radia¬ 
tion is nearly unimpeded (the emerging spectrum is very similar to 
the source spectrum, reported, at peak luminosity, t = 115 Myr, as 
a dashed line) and the continuum normalization drops by ~ 3 (~ 4) 
orders of magnitude in the X-ray (infrared) band. 


3.2 High-redshift MBH mass density 

The approximate flux thresholds for two future-generation observa¬ 
tories (JWST* in the infrared band, ATHENA^ in the X-ray) and for 
the CDF-S^ survey are shown in Fig. 3. We predict that the JWST 
will be able to observe most (> 95%) of the accretion process onto 
a 10^ Mq MBH seed up to a comoving distance corresponding to 
z ~ 25, while ATHENA will only detect ~ 25% of the total evolu¬ 
tion, around the peak luminosity, up to z ~ 15. 

Comparing the predicted peak flux (in the IkeV observed 
band) with the CDF-S sensitivity, we find that this ultra-deep sur¬ 
vey could have observed the accretion process onto a typical MBH 
up to Zmax ~ 15. In the CDF-S survey, Nc = 3 AGN candidates at 
z > 6 = Zmin have been identified (Giallongo et al. 2015, but see also 


* For a NIRcam observation with a Signal-to-Noise ratio of 10 and a total 
integration time of lO'* s. 

^ For a 3cr detection with a total integration time of 3 X 10^ s. 

° Ultra-deep survey in the X-ray with a total integration time of 4 X 10® s. 
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Weigel et al. 2015 where the authors question some of these candi¬ 
dates) inside a sky region of ~ 170arcmin^. An upper limit for the 
number density of MBH seeds can be derived as: 


n.(6 < z < 15) < 


Nc 

ilcDF'Vm- ’ 


(5) 


where SIcdf = 1.1 x 10“® is the sky fraction observed by the CDF- 
S, “y = 2.3 X lO'^Mpc^ is the comoving volume of the Universe 
between Zmin and Zmax^ tD is the duty cycle for accretion and 
T = idethobs is the fraction of time during which the object is de¬ 
tectable within the time frame < z < 15) = 670Myr, assuming 
a single episode of MBH growth at these redshifts^. The predicted 
flux, within the spectral range 0.5-2.0keV, is above the CDF-S 
sensitivity between ~ 85Myr and ~ 115Myr, so that t^et = 30Myr 
and P = 0.045. Using these values we obtain the following upper 
limit for the number density of MBH seeds: 


n.(6 < z < 15) < 


2.5 X 10^5 
D 


(T) 


Mpc 


( 6 ) 


Considering MBHs of initial mass 10^ Mq growing up to ~ 10^ Mq 
(P acucci & Fen'ara 2015), we finally obtain the following upper 
limit for the MBH seeds mass density: 


p.(6<z< 15) < 


2.5 X 10^ 

D 


(T) 


MoMpc 


(7) 


The discussion so far has been limited to the scenario in which 
a MBH with an initial mass close to the peak of the birth mass 
function devised in Ferrara et al. (2014) accretes gas from a LDP 
host halo in the Eddington-limited mode. However, this is not the 
only possible scenario. Table 1 provides a general outline of the ac¬ 
cretion history and CDF-S observability for a MBH seed with ini¬ 
tial mass 10^ Mq in three additional scenarios: standard accretion 
- HDP, slim disk accretion - LDP and slim disk accretion - HDP. 
Table 1 includes the depletion time the duty cycle D, the de¬ 
tection time tdet with its related value of p and the upper limit on 
p,. In the Eddington-limited cases the depletion times are lower 
limits since, while our simulations are run in isolation, in a cosmo¬ 
logical framework the halo growth by mergers and accretion would 
not be negligible within ~ lOOMyr. In the slim disk case, given the 
depletion times of order ~ lOMyr, the hierarchical growth of the 
halo plays a minor role. 

Interestingly, the accretion process is undetectable by the 
CDE-S in the slim disk - LDP case, whose time-evolving spectrum 
is shown in Fig 4. Therefore no constraints on p, can be obtained 
in this case. The slim disk - HDP case is instead observable, albeit 
only for a very short time, fjg, ~ 1 Myr, due to the larger accretion 
rates in this scenario, which produce a sufficiently high luminosity 
despite radiation trapping. The very low value of 'F ~ 1.5 x 10~^ 
produces an upper limit for p, higher than in other accretion sce¬ 
narios by a factor ~ 35. 

Let us now focus on the reason why the slim disk case leads to 
a very different value for p, with respect to the standard accretion 
case. The amount of mass which is available for accretion is equal 
in both cases (~ 10^ Mq) and, to a first-order approximation, also 
the mass actually accreted is similar. This is not strictly true, since 
the radiation pressure is more efficient in the standard case in creat¬ 
ing mass outflows (Pacucci et al. 2015a), but the produced outflows 


'* This calculation assumes spherical symmetry of the host halo and an 
isotropic irradiation of the MBH. If the density along the poles is much 
lower, then a fraction of the sources would be completely unobscured. 


Energy (rest frame) [koV] 

10'^ lO"'' 10'^ 10'^ 10'^ 10° 10* 10* 10° 



Figure 4. As in Fig. 3, but for the slim disk accretion - LDP case. 


are fairly weak regardless of the accretion scenario. The real differ¬ 
ence between the slim disk and the standard cases is due to radiation 
trapping, which decreases the effective bolometric luminosity (i.e. 
the luminosity escaping to infinity) in the former case, with respect 
to the latter. Since the accreting black hole is intrinsically fainter 
in the slim disk case, it will be observable for a smaller fraction 
of time, P (with the extreme P = 0 for the LDP case), during its 
evolution: this eventually leads to a larger upper limit on p,. 

The initial mass of the MBH seed influences the time scale of 
the process and its observability. If, for instance, the high-redshift 
population of MBHs is characterized by a larger average mass, 
the evolutionary time scale would be smaller (since oc oc 
where M^dd is the Eddington accretion rate), but the emitted 
luminosity would be higher (since L oc LEdd 47,). The effect on P 
may vary from case to case, but in general a larger average mass is 
likely to decrease P = tdethobs since tdet ^ tend Mf. This would 
lead to a less stringent upper limit on p,. 


4 DISCUSSION AND CONCLUSIONS 

Using a combination of radiation-hydrodynamic and spectral syn¬ 
thesis codes, we have investigated the time-evolving spectral en¬ 
ergy distribution of an accreting z ~ 10 MBH. The MBH seed, 
whose initial mass is 10^ Mq, is embedded in a dark matter halo 
of total mass IO^Mq. Employing two gas density profiles and two 
accretion modes (Eddington-limited and slim disk) we simulated 
the system until complete gas depletion and we accurately calcu¬ 
lated the time-evolving spectrum of the radiation emerging from 
the host halo. The main results of this work are summarized in the 
following. 

• The spectrum of the emerging radiation, for a MBH observed 
at z = 9, is dominated by the infrared-submm (1 - 1000pm) and 
X-ray (0.1 - 100keV) bands. Photons with frequency shortwards 
than the Lya line are absorbed by the intervening matter at column 
densities Nh > 10^^ cm“^ and reprocessed at lower energies, in the 
infrared band. Due to the very large mean free path Tx > 1 Gpc of 
X-ray photons, growing MBH seeds negligibly contribute to reion¬ 
ization. The continuum normalization is set by: (i) the bolometric 
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Table 1. Accretion history and CDF-S observability for a MBH seed with initial mass M, = 10^ Mq in the four indicated accretion scenarios. 


Accretion scenario 

Obsei'vable 

tend IMyr] 

D 

tdet IMyr] 

T 

p.[MoMpc-3] 

Standard accretion - LDP 

YES 

120 

1.0 

30 

0.045 

<2.5x10^ 

Standard accretion - HDP 

YES 

240 

0.4 

110 

0.16 

< 1.8x 10^ 

Slim disk accretion - LDP 

NO 

7 

1.0 

0 

0 

No constraints 

Slim disk accretion - HDP 

YES 

12 

1.0 

1 

0.0015 

<7.6x10^ 


luminosity of the source and (ii) the column density of the host 
halo. The former determines the overall normalization, while the 
latter determines the ratio between the low-energy and the high- 
energy continua. 

• Our predictions show that the JWST will detect in the infrared 
a fraction > 95% of the accretion process onto a typical MBH seed 
observed at z ~ 9, while ATHENA should observe it in the high- 
energy bands only around the peak luminosity, a fraction ~ 25% of 
the total evolution. Similarly, long-exposure surveys in the X-ray, 
like the CDF-S, could have already observed the accretion process 
on a z ~ 9 object for a comparable fraction of time. The redshift 
of the sources sets their luminosity distance, hence influences their 
detectability. For instance, the standard accretion - FDP system is 
observable by the CDF-S for ~ 41% of the time at Zmin = 6, while 
it becomes undetectable at z ^Zmax =15. 

• From the z > 6 candidates detected in the CDF-S survey (Gi- 
allongo et al. 2015) we estimate the following upper limits on the 
z > 6 MBH mass density: (a) p, <2.5 x lO^MoMpc'^ assuming 
Eddington-limited accretion; (b) p, < 7.6x lO^MoMpc^^ if accre¬ 
tion occurs in the slim disk, highly obscured mode. However, the 
accretion process is undetectable with the CDF-S sensitivity in the 
slim disk accretion - LDP case, due to the flux suppression caused 
by radiation trapping, and no constraints on p, can be given. 

Very recently, we proposed that the first detection of a high- 
redshift MBH seed could have already occurred. Indeed, in Pallot- 
tini et al. (2015) we showed that the observational features of CR7 
(Sobral et al. 2015), a bright Lyu emitter at z = 6.604, may be ex¬ 
plained by accretion onto a MBH of initial mass ~ 10^ Mq. 

For Eddington-limited accretion, our upper limit, p, < 2.5 x 
lO^MoMpc"^, is compatible with the one set by Cowie et al. 
(2012) using observations of faint X-ray sources in the CDF-S, 
while it is more stringent than limits by Willott (2011), Fiore et al. 
(2012) and Treister et al. (2013) (p, < lO^MoMpc^^) and partic¬ 
ularly by Salvaterra et al. (2012) (p, < lO'^MoMpc^^, using the 
unresolved X-ray emission). The current observational constraints, 
however, do not take into account heavily buried, Compton-thick 
objects or radiatively inefficient accretion. Recently, indeed, Co- 
mastri et al. (2015) suggested that recent revisions to the local 
SMBHs mass density, up to p,(z = 0) ~ lO^MoMpc^^, seem to 
imply that a significant fraction of the local SMBHs have grown 
in heavily buried, Compton-thick phases, or by radiatively ineffi¬ 
cient accretion. Our model is a first step towards testing the role of 
Compton-thick or radiatively inefficient phases in the early growth 
ofMBHs. 

In the present paper, we provided a general picture of the inter¬ 
connection between the main accretion mode at work in the high- 
redshift Universe and the black hole mass density. To summarize 
our rationale: assuming that, for MBHs at z > 6, the main accretion 
channel is the standard, Eddington-limited one, we are able to pro¬ 
vide an upper limit on p, which is consistent (and competitive) with 
current estimates. In the standard disk scenario, MBHs have grown 


in heavily buried, Compton-thick phases for about 70% of the total 
evolution time (e.g., ~ 70Myr for the LDP), but they are long-lived 
and intrinsically bright, and visible for a substantial amount of time. 
Assuming instead that at super-critical accretion rates the accretion 
disk thickens, with radiation trapping playing a significant role, the 
suppressed radiative efficiency leads to a much lower intrinsic lu¬ 
minosity. These sources are intrinsically faint, not obscured: the 
Compton-thick phase is short (~ 2Myr, 30% of the total evolution 
time for the LDP) because the obscuring gas is consumed rapidly. 
Our slim disk simulations suggest, as one would expect, that these 
short-lived and fainter MBHs are more difficult to detect in current 
surveys compared to brighter objects accreting in the Eddington- 
limited mode. As a consequence, the upper limit on p, is inevitably 
higher than currently predicted, up to a factor ~ 35. 

MV acknowledges support from a Marie Curie FP7-Reintegration- 
Grant (PCIGlO-GA-2011-303609). GD acknowledges support 
from the Centre Nationale d’Etudes Spatiales (CNES). 


5 SUPPLEMENTARY FIGURES 

In the following page we include additional figures showing the 
time evolution of the emerging spectrum (in v/y units) for the four 
accretion scenarios discussed so far. 
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Figure 5. Standard accretion - LDP case. 
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Figure 6. Standard accretion - HDP case. 
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